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Abstract 

Thermonuclear runaway burning of carbon is in rare cases observed from accreting neutron stars 
as day-long X-ray flares called superbursts. In the few cases where the onset is observed, superbursts 
exhibit a short precursor burst at the start. In each instance, however, the data was of insufficient 
quality for spectral analysis of the precursor. Using data from the propane anti-coincidence detector of 
the PCA instrument on RXTE, we perform the first detailed time resolved spectroscopy of precursors. 
For a superburst from 4U 1820-30 we demonstrate the presence of photospheric radius expansion. We 
find the precursor to be 1.4—2 times more energetic than other short bursts from this source, indicating 
that the burning of accreted helium is insufficient to explain the full precursor. Shock heating would 
be able to account for the lacking energy. We argue that this precursor is a strong indication that the 
superburst starts as a detonation, and that a shock induces the precursor. Furthermore, we employ 
our technique to study the superexpansion phase of the same superburst in greater detail. 
Subject headings: accretion, accretion disks — methods: numerical — nuclear reactions, nucleosyn- 
thesis, abundances — stars: neutron — X-rays: binaries — X-rays: bursts 



1. INTRODUCTION 

Superbursts are rare day-long X-ray flares observed 
from accre ting neutron stars in low-mass X-ray binaries 
(LMX Bs, ICornelisse etaLl [20001: iStrohmaver fc Brown! 
2002]). They are attributed to runaway thermonuclear 
fusion of carbon in an approxima tely 100 m thick layer 
below the neutron star surface (|Cumm ing & B ildstenl 
2001). Superbursts share some properties of the fre- 
quently detected short Type I X-ray bursts (durations 
of ~ 10 — 100 s), which result from runaway burning of 
hydrogen and/or helium accrete d onto the neutron star 
from the binary companion (e.g.. lLewin et al] [l993h For 
example, both have a fast rising light curve, followed by 
a slow decay during which the X-ray spectrum exhibits 
cooling. Superbursts, however, are 1000 times more en- 
ergetic, last 1000 times longer, and o ccur 1000 times less 
frequ ently than the short bursts (e.g.. lKeek fe in 't Zand! 
l2008h . 

Currently 22 (cand idate) superburst s are known from 
13 sources (see, e.g., IKeek et al"1l2012l for an overview). 
All detections are performed with observatories in low 
earth orbit, which implies observations are frequently in- 
terrupted by Earth occultations and passages through 
the South-Atlantic Anomaly. Often the start of a su- 
perburst is not observed. In some cases this makes 
it difficult to prove the superburst nat ure of an event: 
these are referred to as candidates (e.g.. lin 't Zand et al.l 
12001 . Only in 6 cases is there an unambiguous detec- 
tion of the onset. In all these c ases a precursor burst 



is present at the superburst star t (IStrohmaver & Brown 



2002; Str ohmaver fc Markwardtl 1200% lin 't Zand etal 
20031 I2004D . In previous studies the term "precursor 1 



has also been used to describ e a burst that precedes the 
superburst by two minutes ([Kuulkers et al.l l2002bl) , or 
whe n two of these events are vi sible in the superburst 
rise ([Strohmaver fc Brownl [20021 ) . Because these likely 
have a different physical origin, we restrict our study 
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to those precursor bursts that immediately precede and 
transition into the superburst. 

The most detailed precursor observations have been 
performed with the Proportional Counter Array (PCA) 
on the Rossi X-ray Timing Explo rer (RXTE), and are 
from superbursts of 4 U 1820-30 (IStrohmaver fc Brownl 
120021 ) and 4U 1636-53 (| Strohmaver fc Markwardf 2002) . 

4U 1636-53 is a prolific burster that has exhib- 
ited many different modes of t hermonuclear burning 
of hydrogen and helium (e.g.. [Hoffman et al.l 119771 : 
Ivan Paradiis et aLlll986t IRevnivtsev et al.ll2001l). as well 
as carbon burning in 4 s uperbursts ([Wiinandsl 1 2001 : 
I Strohmaver fc Markwardtl l2002t iKuulkers et"aiT 12004 : 
[Kuulkers 200g). 

4U 1820- 3 is a n ultracompact X-ray binary (UCXB, 
IStella et al.l [1987), implying that the material ac- 
creted from the don or star i s likely helium rich 
(Rappapo rt et al.l [T987. see also ICummingl [2003). Its 
persistent flux is observe d to vary on a timescale of 
approximately 171 day s (jPriedhorskv fc Terrell! 119841 : 
iChou fc Grindlavl l200l), and X-ray bu rsts are only de- 
tecte d during periods of low flux (e.g., ICornelisse et al.l 
I2003D. Two (candidate ) supe r bursts have been observe d 
(jStrohmaver fc Brownl 12001 lin 't Zand elaTl I2011bl ). 
The first superburst from this source was observed 
with the PCA, and used to study the changing prop- 
erties of the accretion disk during the superburst 
(|Ballantvne fc Strohmaverl 120041 ) the occurrence of su- 
perex pansion of the photosphere (|in't Zand fc Weinberg! 
2010) , and achromatic variability in the light curve 
in 't Zand et alj|2011al) . 



The precursors resemble short X-ray bursts, and the 
light curve exhibits two peaks, si milar to photospheric 
radius expansion (PRE) bursts (jTawara et al J 11981 
ILewin et al.lll984l ). During powerful X-ray bursts, when 
the luminosity reaches the Eddington limit, the photon 
pressure is sufficient to (temporarily) push out the pho- 
tosphere, causing the surface temperature to drop. The 
lower temperature results in a smaller part of the black- 
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body flux to fall within the energy band in which many 
X-ray detectors are sensitive, producing a dip in the light 
curve. Unfortunately no spectral data of sufficient qual- 
ity have been available to confirm the PRE nature of 
superburst precursors. 

An alternative explanation for the double peaked na- 
ture of precursors has been put forth. Carbon burning in 
a superburst th at proceeds as a deton ation drives a shock 
to the surface ( Weinber g et al.ll2006() . where it produces 
a shock-breakout peak in the light curve, followed by a 
burst from the shock-ignit ed burning of hydrog e n and 
helium in the atmosphere. Weinberg & Bildsten (2007) 
argue that the two peaks are separated by a thermal 
timescale of s econds, similar to the two peaks observed 
in precursors. iKeek fc Hegerl ([201 ID . however, show that 
the fallback of the shock-accelerated atmosphere alone 
dissipates enough energy to power a bright precursor on 
a dynamical timescale of 10 /is. At the same time any 
hydrogen or helium present in the atmosph ere ignites 
and c ontributes to the precursor energetics (IKeek et al.l 
2012). In this scenario the shock breakout and the sub- 
sequent burst are too close in time to be distinguished as 
two peaks, but the burst is bright enough to cause radius 
expansion. 

To test the different models for superburst precursors, 
we devise a technique using anti-coincidence data of the 
PCA to perform time resolved spectroscopy on the two 
superbursts observed with this instrument from, respec- 
tively, 4U 1820-30 and 4U 1636-53. This is the first time 
precursors are studied in detail, and we discuss the im- 
plications for the nature of precursors and of superburst 
ignition. 

2. OBSERVATIONS 

RXTE (|Levine et alJll996j) was launched in December 
1995 and observed t he X-ray sky until January 2012. 
It carries the PCA (jJahoda et al.l [2006) , which has a 
band pass of 2 to 60 keV, and a total collecting area 
of 6500 cm 2 . As part of its science program the PCA 
frequently observed LMXB s, and it detected th ousands 
of short X-ray bursts (e.g.. iGallowav et al.ll2008f ) as well 
as the two mentioned superbursts. 

The PCA consists of 5 identical collimated Propor- 
tional Counter Units (PCUs). Each PCU contains three 
xenon-filled layers, which act as the main counter, and 
are sensitive in the 2-60 keV energy range. On top of 
the xenon detector is a propane-filled anti-coincidence 
detector. It provides additional sensitivity between 1.8- 
3.5 keV with peak sensitivity at 2.5ke\f3, but little cal- 
ibration information is available. Data from the xenon 
layers are typically available per layer and PCU at high 
time resolution and with spectral information, whereas 
for the propane layer only a count rate for all PCUs com- 
bined is available at 0.125 s time resolution without spec- 
tral information. 

At the time of both superburst observations there 
were problems that resulted in a loss of data. At 
the time of the 4U 1820-30 superburst observation, 
RXTE was troubled by a malfunction that prevented 
part of the data from being transmitted to the ground 
([Strohmaver fe Brownll2 002) . During the 4U 1636-53 su- 



perburst the data buffer overflowe d, causing a loss of data 
(|Strohmaver fc MarkwardtJ 120021) . Most notably, high 
time resolution spectral data were lost in both cases. 
During the superburst rises only Standardl and Stan- 
dard2 data are available, which have a time resolution of 
0.125 s and 16 s, respectively. Only the latter mode con- 
tains spectral information, but its time resolution exceeds 
the duration of the precursors. We are, therefore, lim- 
ited to using Standardl data. Additionally, the propane 
count rate is available at 0.125 s time resolution. 

Because of the different energy responses of the xenon 
and propane layers, count rates from the two can be re- 
garded as representing two energy channels. This allows 
for fitting a spectral model with two parameters, such 
as the black-b ody model that is often used to describe 
burst spectra (Sw ank et al.l 119771 ). The p ropane signal, 
howe ver, is not calibrated for this purpose ([Jahoda et al.1 
2006). To investigate the energy response in these two 
"channels" to a black-body spectrum, we perform time 
resolved spectroscopy on several short bursts from the 
two superbursters for which event data with a minimum 
time resolution of 125 [is and with 64 energy channels are 
available. 

We select short bursts from version 0.51 o f the Multi- 
INstr ument Burst Archive^ (MINBAR, IKeek et all 
I2010D . which contains the results of the analysis of 
4,192 bursts observed from 72 sources with the PCA 
([Galloway et al.ll2008T ) and t he Wide-Field Camera 's on 
the BeppoSAX observatory ([Cornelisse et al.l l2003). For 
the purpose of this study we only use PCA bursts. 

The gain settings of the PCUs were changed several 
times during the mission. To be able to compare count 
rates, we restrict ourselves to bursts that were observed 
with the same detector gain as the respective super- 
bursts. Furthermore, the gain was slight ly different for 
each PCU (up to l%. Uahoda et~aT] [20061. The propane 
data, however, are stored as an aggregate, and not per 
PCU. We, therefore, use the propane and xenon data of 
all active PCUs. 

The number and configuration of active PCUs changes 
per observation, and due to the impact of micromete- 
orites PCU and 1 lost their propane lay er on, respec- 
tively , May 13 2000 and December 25 2006 (| Jahoda et all 
2006). We, therefore, report xenon and propane count 
rates per active detector. An active detector that has 
lost its propane layer is still counted as an active detec- 
tor for determining the xenon count rate, but not for the 
propane rate. 

There are small changes in the gain and instrument re- 
sponse over time. For example, xenon is known to have 
leaked into the propane layer. When analyzing the short 
bursts, we generate response matrices for each burst to 
take into account the instrument response and the con- 
figuration of active PCUs at a particular time. 

For spectra and count rates we subtract the back- 
ground measured during a time interval before each (su- 
perburst. This yields the net burst rates and spectra, 
and removes the astrophysical background as well as the 
instrumental background. Here we assume both back- 
grounds remain constant during the burst. We further 
discuss this in Section 14.11 but especially in the burst 



1 We refer to http://heasarc.gsfc.nasa.gov/docs/xte/appendix_f.html 2 See http://users.monash.edu.au/~dgallow/minbar for more 
for XTE Technical Appendix K. details. 
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TABLE 1 

13 PRE Bursts from 4U 1820-30 and a Superburst 
Precursor (Bottom) 



b 5000 
g 4000 



Obs.ID 


Time 


Fluence 


Xenon counts 




(MJD) 


(10 -7 erg 


(10 4 c cm" 2 






cm -2 ) 


PCTJ- 1 ) 


40017-01-24-00 


52794.73813 


3.6 ±0.1 


3.20 ±0.08 


70030-03-04-01 


52802.07557 


4.0 ±0.1 


3.37 ±0.07 


70030-03-05-01 


52805.89566 


4.1 ±0.2 


3.31 ±0.08 


90027-01-03-05 


53277.43856 


3.7±0.1 


3.24 ±0.08 


94090-01-01-02 


54948.82124 


3.2 ±0.2 


3.01 ±0.14 


94090-01-01-05 


54950.70281 


3.5 ±0.1 


3.22 ±0.10 


94090-01-02-03 


54956.77470 


3.3 ±0.1 


3.08 ±0.11 


94090-01-02-02 


54958.73998 


3.2 ±0.2 


2.98 ±0.14 


94090-01-04-00 


54978.32149 


3.7 ±0.2 


3.24 ±0.14 


94090-01-04-01 


54978.49489 


3.7 ±0.2 


3.3 ±0.2 


94090-01-05-00 


54981.18728 


3.8 ±0.1 


3.21 ±0.09 


94090-02-01-00 


54994.53418 


3.2 ±0.2 


2.95 ±0.13 


94090-02-01-00 


54994.61301 


3.6 ±0.1 


3.15 ±0.10 


30057-01-04-08 


51430.07423 


6-8 


4.59 ±0.06 



peak the signal is expected to be dominated by the burst, 
so small changes in the background should not be detri- 
mental to our analysis. 

All extracted data products, including spectra and 
light curves, are corrected for dead time using the pre- 
scription provided by the instrument team. Data re- 
quired for this correction are available at a resolution 
of 0.125 s. When calculating the correction we take into 
account that some of the active PCUs may have lost their 
propane layer. 

3. RESULTS 

3.1. 4U 1820-30 

On 9/9/19 99 the PCA detected a su perburst from 
4U 1820-30 <|Strohmaver fc Brown! [2002]) . The super- 
burst was observed for three hours, and the flux con- 
tinued to decrease when the observation ended. Two 
events with PRE-like peaks and dips are visible in the 
light curve. We consider the first to be the precursor. 
MINBAR lists 14 bursts from this source detected with 
the PCA, 13 of which were observed with the same high 
voltage settings as the superburst (Table [TJ. All bursts 
exhibit PRE. We perform time-resolved spectroscopy on 
the bursts to calibrate the xenon and propane count rates 
such that these can be used to analyze the superburst 
precursor. 

3.1.1. Analysis of short bursts 

Net burst spectra are extracted at a resolution of 
0.125 s in the rise and at the peak of the burst. Af- 
ter the peak we decrease the time resolution by a fac- 
tor 2 every time the count rate per PCU decreases 
by V2j which ensures similar statistics for all spec- 
tra. Each spectrum is fit with a black body in the 
2.5 — 25 keV energy range. Absorption by interstellar hy- 
drogen is taken into account u sing the cross sections of 
iMorrison fc McCammonl (jl983), while fixing the hydro- 
gen column to N H = 0.25 x 10 22 (jAsai et al.ll2000D . The 
fit yields the black-body temperature, kT, and normal- 
ization, K hh = (^) 2 (jn^) , with R the black-body 




lOkpc, 

radius and d the distance to the source. 

For each time bin we extract the count rate in the 
propane layer and in the xenon layer (per PCU). We de- 
termine the sum of the xenon and propane count rate, /, 



10 15 
Time (s) 

Fig. 1 . — For the X-ray burst on MJD 52794 from 4U 1820-30 the 
count rate in the xenon and the propane layers, hardness H (ratio 
of xenon and propane counts), and the results of time resolved 
spectroscopy: the unabsorbed bolometric flux -Fbolj the black-body 
temperature kT, and the black-body normalization K^. For the 
latter we omit two data points with anomalously large values and 
errors. Vertical (green) lines indicate the Icr uncertainties per time 
bin (no uncertainties are indicated for the rates). 

and the hardness ratio of the two rates, H. An example 
of the evolution during a burst of these quantities and 
the spectral fit parameters is shown in Figure [TJ 

4U 1820-30 experiences rather strong radius expan- 
sion during its bursts, which pushes a large part of 
the emission out of the PCA band. The peak of the 
count distribution for a black body lies at an energy of 
£ max ~ 1.59 kT. Because of absorption at low energies, 
the peak of the spectrum is shifted to slightly higher en- 
ergies. When the black body temperature drops below 
kT ~ 1.6 keV the peak is outside the considered band- 
pass. The fits become worse around kT ~ 1.0 keV, re- 
sulting in dips and peaks in the bolometric flux. There- 
fore, our fits are not reliable for kT < 1.0 keV. In further 
analysis, we use kT > 1.0 keV when deriving relations 
between the different parameters. Furthermore, in the 
tail of the bursts the net burst counts are strongly re- 
duced, and the relative uncertainty in H and / become 
large. In our analysis we omit data points where the 
relative uncertainty in H exceeds 50%. 

The hardness ratio traces changes in the spectral 
shape, which for a black body is set by kT. We determine 
the relation between kT and H by fitting a straight line 
(Figure [5]). To determine the goodness of fit we use the 
error in H, because it is typically larger than the error in 
kT. To ensure the latter, we exclude data points where 
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Fig. 2. — Hardness (ratio of xenon and propane counts) versus 
the black-body temperature kT from time-resolved spectroscopy 
of 13 PRE bursts. Crosses indicate the data points for all bursts 
(no uncertainties implied) , whereas disks and diamonds are for two 
bursts on MJD 52794.7 and MJD 54994.5, respectively. The solid 
line is the best linear fit to the data with kT > 1.0 keV. 




(keV) 



Fig. 3. — Black-body temperature kT versus the ratio of the 
black-body normalization i\bb and the total count rate / for 13 
PRE bursts. See Figure [2] for an explanation of the symbols. The 
solid line is the best fit broken power-law to the data with kT > 
1.0 keV. 

the relative uncertainty in kT exceeds 10%. The best fit 
line is 

kT = (0.40 ± 0.06)H keV. (1) 

To obtain \1 — 1 a systematic error of 16% needs to 
be added in quadrature to the uncertainty of the data 
points. This describes the spread in the distribution of 
both individual bursts and all bursts combined. Some in- 
dividual bursts have a smaller spread (Figure [5]), but the 
spread between such bursts is substantial, so this does 
not aid us in improving the fit. Furthermore, the fitted 
line is defined to go through zero. Including an offset, the 
best fit conforms better to the points with kT < 1.0 keV 
(which are excluded from the fit), but the fit statistics are 
not improved; the best fit value is within la consistent 
with 0.0. Furthermore, competition between the fit pa- 
rameters leads to large uncertainties in the parameters. 

Given kT, which sets the spectral shape, is linearly 
proportional to /. We determine the ratio K^/I, and 
we find that as a function of kT it follows a relation that 
is close to a power law (Fig. [3]). We fit a broken power 
law with the break at 2 keV, which is approximately the 
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Fig. 4. — Similar to Figure[T]for the precursor of the 4U 1820-30 
superburst. 

center of the region over which we fit, and we define the 
power-law normalization at the location of the break. 
The best fit is 

^ = (5.0±0.4)xl0-(^)" r , (2) 

with r = 3.8 ± 0.2 for kT < 2kcV and T = 3.10 ± 
0.07 for kT > 2kcV. The power-law index is 1.4 a from 
3.0 for kT > 2keV (where the peak of the black-body 
counts spectrum is within the xenon detector band-pass), 
which is the expected value for a Planck counts spectrum. 
Performing a fit of a single power law with the index 
fixed to this value leads, however, to a larger error in the 
pref actor. 

For all bursts we determine the total number of xenon 
counts per PCU and the unabsorbed bolometric fluence 
(Table [TJ. 

3.1.2. Analysis of the superburst precursor 

We divide the first 13.5 s of the superburst into time 
intervals, using 0.125 s bins at the start to resolve the 
rise and dip in the count rate at the highest available 
resolution. We double the duration of the time inter- 
vals several times to get similar uncertainties in the data 
points of the entire precursor. The tail is resolved at 1 s 
resolution. We extract H and /, and use equations ([!] ^ 
to derive the black body parameters (Figure @| . At the 
start of the precursor we can clearly see increasing, 
accompanied by the dip in kT. This is indicative of PRE. 

The unabsorbed bolometric black-body flux Fboi = 
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Fig. 5. — Similar to Figure Q] for the start of the 4U 1820- 
30 superburst, including the precursor. We omit the data points 
where the xenon rate drops below 0. 

■^-crT 4 can be expressed in terms of I and kT using equa- 
tions (HJ [2]): 

/ kT \ r 

Fbol = / (2keV/ ( 8 - 6±0J ) x 10_12 ergs- 1 cm- 2 , (3) 

with T = 0.2 ± 0.2 for kT < 2keV and T = 0.90 ± 
0.07 for kT > 2 keV. There is a brief dip in F bo i at the 
moment of maximum radius expansion (Figure 0]). This 
is most likely caused by the peak of the Planck spectrum 
being outside the observable band when the temperature 
drops to ~ 1 keV, and not due to a genuine drop in the 
bolometric flux. 

The light curve as shown in Figure [4j i.e. the first 
13.5 s of the precursor, has an unabsorbed bolometric 
fluence of (6.8 ± 1.1) x lO^ergcm^ 2 . After 13.5 s, the 
light curve rises to the superburst peak. If we were to 
consider the precursor as an isolated burst, without the 
interruption by the superburst, the flux would continue 
to decay. We fit F^ i with an exponential, starting at t = 
5.5 s, when the flux d rops below 90% of the peak flux (cf. 
IGallowav~e t al. 2008) . We find an exponential decay time 
scale r = 10 ± 3 s. Using this time scale to extrapolate 
beyond 13.5 s, we estimate an additional contribution to 
the precursor fluence of (1.5±1.0) x 10 _7 ergcm~ 2 , giving 
a total of (8 ± 2) x 10~ 7 erg cm" 2 . 

3.1.3. Analysis of the second PRE phase 

The superburst was only spec trally analyzed from 
100 s after the (precursor) onset bv lStrohmaver fe Brown! 



(|2002T) . lin't Zand fc Weinberg! (|2010l ) study the first 
100 s at a time resolution of 16 s. We use our technique 
to analyze this time interval at higher time resolution 
(Figure [5]). Variability is present in i<boi, most notably 
around t ~ 47 s, which exhibits an anticorrelation be- 
tween kT and ifbb- The variability is present in the 
signal from all active PCUs. 

The superburst shows very strong radius expansion: an 
order of magnitude larger than the radius expansion dur- 
ing the precursor. As suming a distance of 7.6 ± 0.4 kpc 
ijKuulkers et al.l r2003f l . we estimate the velocity of the 
photosphere during expansion by dividing the change in 
radius between consecutive data points by the change in 
time, and find as largest values (2.0 ± 0.7) x 10 2 kms _1 
during the precursor (at t — 0.3 s in Figure [4j and 
(1.6±0.6) x 10 4 kms _1 during the superexpansion phase 
(at t = 14 s in Figure [5]). Note that both velocities are 
each only found from the difference between two data 
points. Velocity measurements directly preceding are 
substantially lower, indicating strong acceleration. Po- 
tentially the velocity of the photosphere increased further 
during the superexpansion, when the signal was lost. 

3.1.4. Comparison of fluence of the precursor and short 

bursts 

IStrohmayer fc Brown! (|2002f ) performed a time- 
resolved spectral analysis of the superburst from 
4U 1820-30 using Standard2 data at 16 s reso- 
lution. They start at 100 s, because at earlier 
times the spectral changes are not well resolved 
by 16 s time bins. At t = 100 s they report 
kT ~ 2.6 keV and K hh ~ 115. We find at that 
time kT = 3.2 ± 0.5 keV and K hh = (7 ± 4) x 10 1 , 
which differ by 1 .1 a and 1.2 a, respec t ively, from the 
measurement by IStrohmayer fc Brown! 120021 (see also 
lin't Zand fc Weinberg! 120101 ). The hardness ratio is 
H = 7.88 ± 0.14, which is at the upper end of the 
interval for which the linear correlation between kT and 
H describes the data well, whereas at higher values of 
H the derived kT values lie somewhat below that line 
(Figure^). 

If we take this as a sign that our derived kT values 
are too large by a factor 3.2/2.6, we can scale equa- 
tion (fTJ) by the inverse of this factor. Redoing the 
analysis with reduced temperatures yields a fluence of 
(5.7±0.8)xl0- 7 ergcm- 2 , and (7.8±1.5) x 10- 7 ergcm- 2 
including an extrapolated tail. The effect on i^bb can be 
estimated from equation ([2]) to be approximately a factor 
2, implying the radius is larger by a factor 1.4. 

The maximum burst fluence we measure is (4.1±0.2) x 
10~ 7 ergcm~ 2 (Table [T]). This is 2.5 a away from our 
original fluence measurement and 2.0 a from the one 
with the scaled kT; here we did not include the esti- 
mate of the extrapolated tail fluence, as this increases 
the uncertainty. Whereas the inherent uncertainty in 
the method we employ does not allow for a strongly sig- 
nificant difference in the precursor fluence and the maxi- 
mum burst fluence, a direct comparison of the integrated 
burst counts clearly shows the precursor to be more en- 
ergetic (Figures [IJ [4j . The normal bursts have at most 
(3.37±0.07)xl0 4 ccm~ 2 (Table[l]), whereas the precursor 
(without extrapolated tail) has (4.59±0.06) x 10 4 ccm~ 2 , 
which is a 17 a difference. The precursor has at least 



Keek 



L 3000 
O 2500 
&h 2000 
7 1500 
"> 1000 
-H. 500 

14 
12 
10 



> 

o 



^ I Propane ^, 



— hju. j^M— — i 




Fig. 6. — Similar to Figure[T]for the precursor of the 4U 1636- 
53 superburst. Data points with particularly large uncertainties in 
the hardness are omitted. 



1.36 ± 0.03 times more counts, and between 1.4 ± 0.2 
(scaled, no tail) and 2.0 ± 0.5 (original, with tail) more 
fluence. From the difference in count rates, we conclude 
that the precursor is significantly and substantially more 
energetic than the normal bursts. 

3.2. 4U 1636-536 

On 2/22/2001 the PCA observed a superburst from 
4U 1636-53. The light curve exhibits a double peak at 
the onset (Figure [5]). The number of active PCUs with 
a propane layer was the same as during the superburst 
from 4U 1820-30 (3), but the observed propane count 
rate was less than half. The uncertainty in the hardness 
ratio, H, is large in the first few seconds. Coincident 
with the dip in the rate, a hint of a dip in H may be 
present around 1 s after the onset, indicating a dip in 
the black-body temperature, which would be consistent 
with PRE. The uncertainties, however, are far too large 
to draw any conclusions. We applied the same analysis 
as for 4U 1820-30, using three PRE bursts detected in 
the preceding month with the PCA (Table |5J). We find 
relations analogous to Equations (Q] [2]) : 



kT = (0.307 ± 0.012)iT keV, 



(4) 



which requires a systematic uncertainty of 7% to obtain 
xl = 1, and 



(4.9 ±0.7) x 10" 



kT V 
2keV/ 



TABLE 2 

3 PRE Bursts from 4U 1636-53 and a Superburst 
Precursor (Bottom) 



Obs.ID 


Time 


Fluence 


Xenon counts 




(MJD) 


(10 -7 erg 


(10 4 c cm" 2 






cm -2 ) 


PCU" 1 ) 


50030-02-04-00 


51937.11615 


4.93 ±0.12 


4.25 ±0.09 


50030-02-05-00 


51942.10027 


4.78 ±0.14 


4.49 ±0.10 


50030-02-05-01 


51941.87561 


4.39 ±0.10 


3.87 ±0.08 


50030-02-08-01 


51962.70299 


3.5 ±0.9 


2.69 ±0.06 



(5) 



with T = 4.2 ± 0.3 for kT < 2keV and T = 2.6 ± 0.3 
for kT > 2keV. These values are consistent with those 
obtained for 4U 1820-30. Using these relations we deter- 
mine the black-body parameters for the superburst pre- 
cursor. There is a hint of weak radius expansion, but the 
errors are too large to constrain whether PRE is present 
(Figure [6|) . 

Because the rise of the superburst emission in the light 
curve is slow compared to the superburst from 4U 1820- 
30, it is more difficult to separate the precursor flux from 
the superburst flux. We determine the total counts and 
the fluence up to the local minimum in the flux at t — 
18 s (Table [2]). For the three short bursts, the counts 
and fluence at t > 18 s contribute at most 15% of the 
total. Furthermore, part of the counts and fluence at 
t < 18 s may be attributed to the superburst. Therefore, 
we regard the counts and fluence up to t = 18 s, increased 
by 15%, to be an upper limit to the precursor emission. 
Because of the large uncertainty in the flux, the fluence 
is not significantly different from the fluence of the short 
bursts, but the counts are significantly lower. 

4. DISCUSSION 

We devised a method to perform time-resolved spec- 
tral analysis of superburst precursors using data from 
the anti-coincidence detector of the PCA on RXTE. For 
a superburst of 4U 1820-30 we clearly see PRE, but 
for 4U 1636-53 the data quality is insufficient to argue 
whether radius expansion occurred. In this section we 
first discuss the possibility for improvements to our spec- 
tral analysis technique. Then we discuss the implications 
of our results. 

4.1. Improvements to the analysis method 

The spect rum of Type I burs ts is typically well fit by a 
black body (|Swank et al.lll977l ). Using the propane and 
xenon counts we have precisely enough energy "channels" 
to fit this model. There are improvements over a pure 
black body that we cannot implement. Corrections are 
proposed in the form of a so-called color correction factor 
for the black-bo dy temperature of typica lly several tens 
of percents fe.g.. ISuleimanov et al.ll201lD . Furthermore, 
we obtain the net burst spectrum by subtracting the per- 
sistent background from a pre-burst time interval. Here 
we assume that the pre -burst emission is dom inated by 
the accretion disk fe.g.. iKuulkers et~a l. 2002a), and that 
this emission remains unchanged during the burst. When 
there is strong radius expansion this may no longer be 
correct. Indeed, the count rate drops below the pre-burst 
level during the long PRE phase in the superbu r st rise 
of 4U 1820-30 (Figure©. IStrohmaver fc Brownl ((2001) 
find that in the tail the qualitative results do not differ 
greatly between assuming a constant persistent flux and 
including the persistent spectrum in the fit. 
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iStrohmaver fc Brownl (|2002f) include in their spec- 
tral model for the superburst of 4U 1820-30 an emis- 
sion line and an absorption edge. These features pos- 
sibly result from reflection of burst emission off the 
accretion disk. Later studies have successfully used 
spectral models including a dis k reflection component 
(iBallantvne fc Strohmaverl 12004 lin't Zand fc Weinberg! 
120101 ). We do not include these in the model when fitting 
the normal bursts. Even though they reach the Edding- 
ton limit, the radius expansion is smaller for the normal 
bursts than during the superburst. Furthermore, during 
the superburst longer time intervals can be chosen for the 
spectral analysis, which allows these spectral features to 
become significant. This is not the case during the short 
bursts. 

There is considerable spread in the relation between 
the hardness ratio and the temperature (Figure [2]). The 
width of the distribution is due to the spread in the 
data points of individual bursts and the spread between 
bursts. We investigate the possibility of improving on 
the fitted linear relation. The xenon detector's effective 
area is well known, but the area is poorly constrained for 
the propane layer. An effective area curve, that has the 
shape of a power law with a smooth low-energy cut-off, 
is presented in the XTE Technical Appendix F. We em- 
ploy this curve and a black-body model in an attempt 
to reproduce the propane counts for the short bursts. 
Whereas at low count rates there is reasonable agree- 
ment, especially at the higher count rates during the 
burst peaks there are strong deviations from the observed 
count rates, with a large spread in the deviations. The 
effective area curve is, therefore, unusable to improve the 
H — kT relation. 

4.2. Precursor from 4 U 1820-30 

The start of the PCA superburst from 4U 1820-30 ex- 
hibits two instances of a pea k and dip in the light curve 
(|Strohmaver fe Brownl I2002T) . The first we refer to as 
the precursor, whereas the second has been identified as 
a period of superexpansion (jin't Zand fe W einberg 120101 
see Section |4~3|) . 

The precursor clearly exhibits PRE (Figure |4} . The 
peak of the black-body temperature is reached 4 s after 
the onset (the "touch-down point"). The peak values of 
the black-body parameters are similar to those of the 
short bursts (Figure [lj. The precursor has, however, a 
longer decay time and a larger fluence. It is more ener- 
getic by a factor 1.4 — 2 than the most energetic of the 
short bursts we analyzed (Table [T]). 

If a sufficiently large column of carbon ignites, burning 
is pre dicted to proceed as a detonation ([Weinber g et al.l 
I2006I). A shock is generated, and travels to the sur- 
face ([Weinberg fc Bildstenl 12007). This produces a brief 
shock breakout peak in the light curve. The shock 
heats the atmosphere, which produces a precursor burst. 
Furthermore, it triggers the burning of hydrogen and 
heli um, which adds to the en er getics of the p recur- 
sor ([Weinberg fc Bildsted[200l . feeek fc Hegen (pOll 
showed that most of the heat from the shock is re- 
leased upon the fallback of shock-accelerated layers. 
At that moment any available hydrogen/helium ignites 
(|Keek et al.ll20 12h . The heat released by the shock is ex- 
pected to be of similar magnitude as that produced by 
helium burning in the atmosphere. 



The superburst took place at a level of persistent 
flux where short bur sts are observed from this source 
(jZdziarski et al.lfeOOTf ) . The amount of helium present in 
the neutron star atmosphere at the moment of the super- 
burst depends on the time since the last occurrence of a 
short burst. Unfortunately this burst was not observed, 
but the fuel column is at most of similar size as the most 
energetic short burst. 

Given that we observe the precursor to have a larger 
fluence than the short bursts, the burning of helium is 
insufficient to power the precursor. We propose shock 
heating as the most likely source to complete the pre- 
cursor energetics. This provides a strong indication that 
the superburst indeed proceeded as a detonation, and 
confirms the scenario of a shock-induced precursor. 

The two peaks in the 2 — 60 keV precursor light 
curve are caused by PRE; they are not two peaks 
intrinsic to the (bolometric) signal. This disproves 
the suggestion that the first peak is caused by shock 
breakout and the second is powered by a helium flash 
(jWeinberg fc Bildstenl [2001 . A scenario where both 
these events are separated in time by an interval that 
is shorter than the 0.125 s resolution of the data, such as 
the dynamical ti mescale of ~ 10~ 5 s, is consistent with 
the observations (|Keek fc Hegerll2011h . 

4.3. Superexpansion phase 

The 4U 1820-30 superburst is one of only few bursts 
to exhibit superexpansion, where briefly the black- 
body radius increases by a s much as a factor 100 
(|in't Zand fc Weinberg! 120101 ). It is speculated to be 
caused by the ejection of a shell of material from 
the neutron star surface. The superexpansion disap- 
pears from the observations when the shell turns op- 
tically thin, and a lower l ying layer becomes visible 
(jin't Zand fc Weinberg 2010). This lower layer is moder- 
ately expanded and is observed to settle over time. This 
superburst has the longest period of moderate expansion 
observed to date of close to 1400 s. 

Using our method we can study the superexpansion 
of the superburst at much higher time resolution (Fig- 
ure Et . Using a distance of 7.6 kpc (jKuulkers et alJ 
2003), we find black-body radii of up to 280 km, and 
we measure the rate at which the black-body radius in- 
creases to be (1.6 ± 0.6) x 10 kms -1 , which are similar 
to the values found for two o ther superexpansion bursts 
(jin't Zand fc Weinbergl l2010). Including corrections can 
lead to 40% higher values (Section 13.1 4j) . The rate of 
black-body radius increase may trace the velocity of the 
expanding photosphere, but deviation of the spectrum 
from a black body and the changing optical depth have 
to be taken into account to determine the true photo- 
spheric velocity. Furthermore, during the strong radius 
expansion phase variations in the persistent flux may be 
present, which have not been accounted for. 

4.4. Variability in the superburst light curve 

Achromatic variability is present in the tail of the su- 
perbur st of 4U 1820-30 between 3 365 - 7200 s after the 
onset (jStrohmayer fc Brownl 120021 ) . This variability is 
suggested to result from alternate blocking and reflect- 
ing of burst emission by density inhomogeneities in the 
accreti on disk, that are induc ed in the superexpansion 
phase (jin 't Zand et alj|2011ah . 
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Shortly after the superexpansion period (around 47 s 
after onset or 33 s after the onset of superexpansion), an- 
other peak is present in the flux. The peak is also visible 
in the hardness ratio, so this is not achromatic variability 
and must be of a different origin. We also find the vari- 
ability in the black-body parameters. These were, how- 
ever, derived assuming a black-body spectrum, whereas 
the variability need not be of this nature. For example, 
the persistent background may vary, whereas we a ssume 
a constant background, lin 't Zand et al.l (|2011al ) esti- 
mate different timescales for processes in the disk. The 
timescale closest to that of the observed variability is the 
viscous timescale of approximately 60 — 950 s on which 
the accretion disk settles. 

4.5. Precursor from 4 U 1636-53 

The count rate for the superburst precursor of 
4U 1636-53 exhibits a dip similar to the light curve of the 
4U 1820-30 superburst. Only a weak hint of PRE, how- 
ever, may be visible in the hardness ratio. The quality 
of the data at the onset of the superburst is insufficient 
to argue for the presence of absence of radius expansion. 
The total photon counts are significantly lower for the 
precursor than for three short bursts in the preceeding 
month. Unlike the precursor of the 4U 1820-30 super- 
burst, therefore, this precursor is not more energetic than 
what c an be provided by nuc lear burning in the atmo- 
sphere. ICummirig~e t al. (2006j) derive for this superburst 
an ignition depth that is j ust above the minim um re- 
quired to generate a s hock (Keck & Hegcr 2 Oi l; see also 
iWeinberg et al.ll2006ft . and the models of iKeek fc Hegerl 
(|2011h that best fit the light curve do not produce a shock 
at all. The precursor would, therefore, mostly be pow- 
ered by the burning of the available hydrogen and helium 
in the envelope, the amount of which depends on the time 
since the previous short burst, which is unknown. This 
may explain the lack of a powerful precursor and a strong 
PRE phase. 

4.6. Outlook 

The time resolution of 0.125 s was insufficient for fully 
resolving the start of the radius-expansion phase of the 
precursor from 4U 1820-30, and the collecting area was 
too small to study the 4U 1636-53 precursor. A future 
mission such as the Large Observatory for X-ray Timing 



(LOFT. iFeroci et al.ll20Tll) is required to study precur- 
sors in detail. This will allow us to confirm the PRE 
nature in other superbursts, and provide important in- 
formation on the detonation process deep in the neutron 
star envelope. 

5. CONCLUSIONS 

Using data from the propane detector of the PCA on 
RXTE, we perform the first detailed spectral analysis 
of superburst precursors. For the PCA superburst from 
4U 1636-53 the data are of insufficient quality, but for the 
PCA superburst from 4U 1820-30 the precursor clearly 
exhibits PRE. This confirms that the precursor is a single 
peak in the (bolometric) emission, and the dip in the 
peak is due to the limited band pass of the instrument. 

The precursor from 4U 1820-30 is 1.4—2 times more 
energetic than the most powerful short X-ray burst we 
analyzed from the same source. This shows that the 
thermonuclear burning of helium in the atmosphere is 
insufficient to power the precursor. We suggest shock 
heating as the most likely additional contributor to the 
energetics. This is strong support for recent numerical 
models that predict superbursts to proceed as a detona- 
tion, and that the generated shock deposits enough heat 
in the atmosphere to power a bright precursor burst. 

With the same technique we study the superexpansion 
phase of the 4U 1820-30 superburst in greater detail then 
was previously possible. We find similar expansion fac- 
tors as reported for the few other superexpansion bursts, 
and we derive expansion velocities of up to 5% of the 
speed of light. 
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International Team on Type I X-ray bursts. The author 
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